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We address the question of whether an anisotropic gap d^2_y2 symmetry is compatible with 
locahzed states in the normal phase. The issue is important in high Tc superconductors where a 
superconductor to insulator transition is observed in the underdoped regime, together with a number 
of experiments that support d-wave pairing. We find a reentrant behavior of superconductivity in 
the strongly disorder phase. 
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. There is a growing body of experimental evidence in high Tc superconductors that indicate that the pairing state is 
0^ ■ of d^2_y2 symmetry [Q. In superconductors with an anisotropic order parameter, both magnetic and non-magnetic 
impurities are pair breaking. For d-wave symmetry, the effect of non-magnetic impurities is equivalent to magnetic 
impurities in s-wave superconductors |^ . Perturbation theory in the impurity scattering introduces a mean free path 
C I £ for the extended states, and the standard treatment indicates that anisotropic superconductivity is destroyed when 
~ 0.1 1^, with ^ the coherence length. The situation is somewhat different for the so called extended s-wave 
symmetry. This corresponds to an order parameter with uniform sign which could, in particular, vanish at certain k 
directions B . Point impurities are not pair-breaking in this case, but they are "pair-weakening" : for small impurity 

L . . . .... . . (cr) 

^-H , concentration nimp, 7c decreases linearly with riimp, but the critical impurity concentration is nl^^ oo. On the 

^ 1 other hand, charge dynamics in oxide superconductors is basically two-dimensional, and it is known from the scaling 

C [ theory of localization that in two dimensions all the one particle states are localized Q . Actually, the experimental 

. evidence from resistivity measurements for low-doping is consistent with a divergent resistivity as T — > that is 

' cut-off by the superconducting transition [|] . The resistivity shows an upturn at a characteristic temperature T^in 

, that is apparent when T^, < T^i^. Qualitatively, T^i^ corresponds to the temperature scale for which the inelastic 
scattering length is comparable to the localization length. Conversly, if the elastic mean free path is much bigger than 

^\ the coherence length, then Tc > Tmin the localization effects are not important, and the variation of Tc with disorder 

"Y^ ■ will be given by the usual pair breaking expressions. In the regime of Tc < Tmin, it is clear that superconductivity 

Ci ' establishes at a temperature low enough for the effects of localization to be of importance. The purpose of the present 

^ i work is to present a treatment of anisotropic superconductivity that incorporates the fact that states from which 

I , the superconducting state is built up are localized, and reconcile two seemingly conflicting properties: the observed 

^ , insulator-superconductor transition, and anisotropic pairing. We also discuss a new prediction that emerges from our 

g theory regarding reentrant superconducting to insulating transition at low temperatures. 

^ , For concreteness we consider fermions on a lattice described by the following Hamiltonian 

.si: H = H,-UY,K^^, (0.1) 

' with Hq being a on e-electron Hamiltonian that includes disorder, with eigenstates ipu{x) and eigenvalues e^. The 



second term in (0.1) corresponds to an instantaneous attractive interaction with an implicit cutoff at a characteristic 
energy lud. In order to model d^2_y2 symmetry we choose Aj, of the form 

= 7f E {^lA+si - <^ld+s^) ' (0-2) 

with 6 = ±ei,±e2 being the lattice versors, and e±ei = — e±e2 = 1- We argue below that the effects of localization 
on the critical temperature Tc for d-wave pairing are qualitatively the same as those for p-wave pairing, an example 
of which is the spin polarized phase 0, that we model with AJ,^ = ^s'^xa-'^t+Sa' ^+ei ^ ^^-ei = I7 and e+e2 — 
e-ea = 0. 

The critical temperature is determined by the self consistent solution of the following linearized gap equation Q 

x,x')A,,>, (0.3) 

X 

where A^^ = (A^;), and the kernel K, written in terms of the exact eigenstates of Hq is given by 
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^(-. - ) = 2 E E ^^^^ (e.-z.„)(e,+.c.„) ^^'^^ 

with T the temperature and uj^ — (2rt + l)7rT the Matsubara frequencies. Also, we have taken h = ks = I- From 
now on we will take the eigenstates as real. 

We next assume that the gap is uniform, A^,- = |A|, which is justified for ojd ^ AT4^, with AW the typical level 
spacing between states within a localization length of each other. In that case we can integrate (^^) over x and x' 
and reach the condition 

with Np the density of states at the Fermi level, and 5(1^) is the power spectrum of the operator D = 

e<5 {\x + S){x\ + \x) {x + 6\) ,given by 



giu) = y\{f,\D\,y)\^d{s. - SF - col (0-6) 



where the symbol represents an average over states /i such that — ep- For s-wave symmetry, D corresponds to 
the density operator D — |a;)(a;| , and = since the density response is not sensitive to scattering (in the 
g = limit, which is our case of intere st). This is valid even when the states tpu{x) are localized. With this frequency 
dependence of g{cj) replaced in Eq. ( |0.5| ) one obtains an equation for the critical temperature identical to that of 
the pure system. This is the extension of the Anderson theorem to the case of localized states, which was discussed 
by Ma and Lee using a variational approach We conclude that under the above assumptions {lod ^ AW and 
uniform gap), the effects of localization on the critical temperature are contained, through Eq.( |0.5| ), on the frequency 
dependence of spectral function of the operator that has the symmetry of the order parameter. The function g{ijj) 
can be calculated diagramatically, since it is given by a two-particle bubble with bare vertices 7^ = cos kx — cos ky 
for d-wave, and 7^ = isinfe^, for p-wave symmetry. From now on we ignore the lattice effects, and take 7^ = cos20kj 
which corresponds to a ga p function A(k) = A(T) cos 20k- 



We first write Eq. (0.6) as 



giuj) = Re V cos 29 (^) cos 20k' , (0.7) 

with 

$k,k' = (G« (k, k'; £f )GA(k', k, sp + uj)) , (0.8) 

where now (. . .) denotes impurity average. We follow the work of VoUhardt and Wolfle (VW) ML and compute g{uj) 
within the self-consistent theory of localization. We prove that the frequency dependence of ( |0.7| ) is essentially the 
same as that of the conductivity. The quantity $k,k' obeys the Bethe-Salpeter equation 

$k,k' (w) = G^G^ <5k,k' + G£G^ ^k,k" (c^)$k",k' (w) , (0.9) 

k" 

with C/k,k"('-^) the irreducible vertex function. In Anderson localization, smi^/e-particle quantities (e.g. the density 
of states) are smoothly varying functions of disorder. It is then reasonable to approximate the self energy E by the 
lowest order result in the impurity scattering Uq, and use the Green's functions in the form 

G^{e) = , (0.10) 

e - ek + 27 

with, 1/t = 2ttNpUo, and G^ — (G^)*. We have assumed a S correlated disordered potential u{x), such that 
{u{0)u{x)) — UQn\n^p6{x) = UoS{x). Using this expression for the Green's functions we can write the Bethe Salpeter 
equation as a kinetic equation in the form 



(w - -)4'k.k' = -AGk 

T 



5k,k' +E^k,k"«'k" 



(0.11) 
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with AGk = — G^. If we replace in (3.11) the irreducible vertex by the bare vertex Uq, we obtain 



A very simmilar expression was obtained before in a treatment of the Raman response in the I = 2 channel in 
the presence of impurities However, since the Ram an re spon se R jtj) is given by a correlation function, there is 
an additional factor of uj and RUli) — Log{uj). Inserting ( |0.12| ) in ( |3.5| ), we obtain the well known expression for the 
critical temperature variation - ln(T'c/Tco) = ^'(1/2 + 1/AtttTc) - -0(1/2)- 

Following VW, we observe that since AGk is srongly peaked a,t k = kp, the dependence of <I'k,k' on the magnitude of 
the wave vectors will be dominated by AGk- We define $k — X^k' '^k,k' cos26'k', and extract the angular dependence 
in k using a Legendre expansion in which we keep up to the I — 2 term: 

AG 

$k = — — ^ V [1 + 2 cos 6lkCos6'k' + 2cos26lkcos26'k'] «'k' - (0-13) 

— ZTTiNp ^-^ 
k' 



Multiplying ( D.ll ) by cos26'kCOs26'k', summing over k, k', and using ( |0.13 ), we obtain 



V cos 20k$k,k' cos 20k' = ^^^^ , (0-14) 

^ uj- M(uj) 

where M{u!) is a '7 = 2 relaxation kernel", given by 

M{uj) = - + -4r- V cos 26'kAGkC/k k' AGk' cos 26ik' - (0-15) 

T TTjyp ^ — ' 
k,k' 

The structure of M(lu) is very simmilar to that of the current relaxation kernel- Note that in deriving this expression 
we have neglected terms that mix different angular dependences in M{uj), and that give rise to factors cos mO]^ cos 26'k' , 
with m = 0, 1, 2. These terms do not appear in the treatment of VW. Neglecting these terms is justified, since we are 
anticipating the inclusion of the contribution to the irreducible vertex that originate the divergence of M(lu). The 
infrared divergence of M comes from the maximally crossed (MC) diagrams ||l^, which are irreducible, and contribute 
with 

for k ~ — k', and with Dq the bare difussion constant. Due to this divergence we can take k = k' for t he angular 
integral and the "off-diagonal" contributions vanish due to orthogonality. Using the above U^^, in ( |3.15 ) we obtain 



the logaritmic low-frequency divergence M{uj) ^ logu) that is familiar from perturbation theory of the conductivity 
[0. Since the low frequency I = 2 kernel is essentially the same as the current relaxation kernel (the angular 
integrations give the same result), it will still be related to the difussion constant D{u!) ~ iDo[AI{uj)T]^^ . This allows 
us to go beyond perturbation theory and determine M (lu) self consistently through the equation 

T T UJ — k'^UQlAl (uj)t\ ^ 
Equation (|0.17D can be solved for low frequencies, giving 

M{uj) = 1 - (0.18) 



and 

47r 1 + (a;rj^ 

with UJ ~ UJ + ujq/uj. The characteristic frequency c^o is finite in the localized phase, and is given by ujq — vp/{\/2X), 
with lambda the localization length. This result implies an expression for the dependence of the critical temperature. 



3 



which can be considered as a generahzation of the Abrikosov-Gorkov-Maki (AGM) formula for locaUzed anisotropic 
superconductors. Our result is then 



ln( — ) = -ln( )— + — 

Ic UJD T T+ 



p+) 



T 



(0.20) 



where l/r± = ^/w^ + (1/2t)2 ± 1/r, and p± = l/47rr±Tc. 

The above expression gives the change in critical temperature as a function of t and ojq- In contrast with the AGD 
formula, in our case the relative change in Tc is dependent of the cutoff frequency ujd- In order to visualize the effects 
of localization on Tc predicted by our treatment, in Figure 1 we show plots of Tc as a function of r, treating lvq as 
an independent parameter. In principle these two quantities are not independent, but the plot emphasizes the fact 
that when Tc/Tco < 0.1 there is a reentrant behaviour to a normal phase at very low temperatures. We note that 
the reentrant behaviour correponds to a parameter range that is realistic in high Tc oxides. The curves of Figure 1 
show that there is reentrant behaviour when ujqt ~ 1. From the self consistent theory, ujot — ^/\/2A, where A is the 
localization length. On the other hand j = (e^'^^^ — 1)^/^. In high Tc oxides, since the particle density is small, kpi 
can be of order one, and u>ot ~ 1. We conclude that localization effects are important in high Tc oxides, and that 
d-wave superconductivity can coexist with Anderson localization. 

Finally, we discuss briefly the mechanism underlying the reentrant behavior in Tc. The reentrance to a normal 
state as temperature decreases can occur if the entropy of the superconducting phase is higher than the entropy of 
the normal state at very low temperatures. This can happen if the localization length is of the order of the mean free 
path. In this regime one can estimate the density of states p{E) of the quasiparticles by computing the quasiparticle 
energies Ei, as a correction of the particle energies in second order in the pairing interaction jsj. Since states that 
are degenerate have negligible overlap in the localized phase, we can extract the typical correction to the quasiparticle 
energy from (/(w), and obtain 



piE) 



No 

I dE I 
I de I 



No 



1 



de J e+u 



(0.21) 



The result is then that the density of states at the fermi energy can be bigger than the normal state density of states, 
and the superconducting entropy is higher at low temperatures. The behavior of p{E) is nonmonotonic in E and on 
e gets two transitions as temperature increases. 

In comparing with experiments one should be able to isolate the effects of hole doping and the effects of disorder. 
Recent experiments on Bi2Sr2Ca(Cu,Co)08+j^ [|l5| it is shown that Tc shows a drastic decrease with increasing disorder 
and the resistivity shows a minimum suggesting Anderson localization. From angle resolved photoemission it is 
concluded that the carrier concentration is unchanged, as expected since Co is a 2"*" ion. Hall effect measurements 
are also indicate that the carrier concentration does not change ||l^. From Ma and Lee's theory ||] one does not 
expect a big decrease in Tc if the symmetry is s-wave, since Anderson theorem still holds if the mobility edge 
crosses the chemical potential. The fact that Tc is changing in the "coexistence region" between the insulating and 
superconductivity phases, can be better understood if the symmetry is d-wave. The reentrant behaviour has not been 
measured to our knowledge, and it constitutes a well defined experimental test of the present theory. 

In summary, we have shown that d-wave superconductivity is compatible with Anderson localization of the one- 
particle states. In addition, out treatment extends previous calculation of the Raman response in disordered systems 
to the localized phase. We also made the interesting prediction of reentrance to the normal phase as temperature is 
lowered in the superconducting phase, when disorder is strong. 
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FIG. 1. Critical temperature vs. inverse scattering time F = 1/4-ktTco for fixed values of wo = 0. (solid), 0.1 (dashed), 0.2 
(long dashed) 0.25 (dot-dashed). 
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